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Abstract 

It is known that hermitean random matrix ensembles can be identified with symmetric 
coset spaces of Lie groups, or else with tangent spaces of the same. This results in a 
classification of random matrix ensembles as well as applications in practical calculations 
of physical observables. In this paper we show that a large number of non-hermitean 
random matrix ensembles defined by physically motivated symmetries - chiral symmetry, 
time reversal invariance, space rotation invariance, particle-hole symmetry, or different 
reality conditions - can likewise be identified with symmetric spaces. We give explicit 
representations of the random matrix ensembles identified with lateral algebra subspaces, 
and of the corresponding symmetric subalgebras spanning the group of invariance. Among 
the ensembles listed we identify as special cases all the hermitean ensembles identified 
with Cartan classes of symmetric spaces and the three Ginibre ensembles with complex 
eigenvalues. 



1 Introduction 



Ensembles of random matrices appear with increasing frequency in various branches of 
physics. In a typical application, an appropriately chosen ensemble of random matrices 
(hermitean or not) replaces an ensemble of physical quantum operators with known global 
symmetries in problems involving highly complex, interacting systems. The method, pio- 
neered by Wigner |1] and Dyson [2j over 50 years ago in a statistical approach to spectra 
of nuclear energy levels, can produce useful information of a universal character on the 
spectrum of eigenvalues of such operators, ranging from Hamiltonians to transfer matri- 
ces, scattering matrices, and Dirac operators. The most characteristic feature of spectra 
described by random matrix theory is the presence of strong geometrical correlations be- 
tween eigenvalues. Dyson p], followed by Hiiffmann [4], were the first to point out the 
connection between ensembles of random matrices, subject to certain symmetry condi- 
tions, and geometrical manifolds referred to as symmetric spaces. Further work on this 
was done by Altland and Zirnbauer and by Caselle in references O [6l [7]. In reference [8] 
the classification of symmetric superspaces in random matrix theory is discussed. 

The relative simplicity with which universal spectral characteristics in widely diverse phys- 
ical systems can be described in this approach, has led to a huge increase in activity in 
the past two decades or so, and the literature in the field by now comprises thousands and 
thousands of pages (for reviews see [9]). Initially only hermitean operators were studied, 
and ensembles were classified according to Cartan's scheme of symmetric spaces, using a 
group theoretical viewpoint (for a review of this approach see [IQJ). Because of the iden- 
tification of the set of random matrix eigenvalues with radial coordinates on a symmetric 
space, the distinguishing characteristics of the matrix ensembles follow naturally via the 
Jacobian from the underlying geometry of the coset space and the root system of the asso- 
ciated Lie algebra subspace. This geometrical connection has so far been applied mainly in 
condensed matter physics [HI [HI 112] , but the possibility of applications in other branches 
of physics should not be overlooked. 

We observe that most of the time, the diffusion operator on the symmetric space, cor- 
responding to the first Casimir invariant associated to the enveloping algebra, played an 
important role in these applications. This was the case also in |T2], where physical mea- 
surable quantities of a carbon nanotube were expressed entirely in terms of quantities 
characterizing the root lattice of the sjTiimetric space on which the transfer matrix lives 
(for details see section [2]). This shows that the connection with an underlying symmetric 
space is not just a theoretical concept, but can be used for determining physical observables. 

Lately however, an increasing number of authors employ non-hermitean random matrix 
ensembles in a range of different contexts (for a review with emphasis on applications in 
quantum chaotic scattering, see |13]). Except for being interesting in their own right. 
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these new ensembles are relevant to unexplored phenomena in field theory and condensed 
matter physics, and appear also in other contexts, for example in the dynamics of neural 
networks. Whereas the range of eigenvalues of a hermitean ensemble is found on the 
real axis, the corresponding support of eigenvalues of a non-hermitean matrix ensemble 
occupies a two-dimensional domain in the complex plane. In this paper we take a first step 
towards clarifying the issue of whether there is a corresponding geometrical structure that 
characterizes the relatively less known non-hermitean "relatives" of the known, classified 
hermitean ensembles, taking as a starting point the paper by Bernard and LeClair [H] 
in which non-hermitean random matrix ensembles were classified according to symmetry. 
Except for two new ensembles apparently forgotten in [H], we verify the classification of 
this reference. The main result of this work is to show that all the resulting ensembles 
correspond to geometrical symmetric spaces, the classification of which goes beyond the 
standard Cartan classification of hermitean random matrices (we do expect, however, that 
the non-hermitean ensembles will be identified with Cartan classes possessing a complex 
structure). As we explain later, we do not claim that the list of symmetric spaces is 
complete, nor that it comprises all possible random matrix ensembles. However, it does 
appear to include many of the purely random matrix ensembles studied in the literature. 
The most important exception arises if one adds a fixed matrix to a random ensemble, or 
imposes an extra condition like the condition of normality (a matrix N is called normal 
if [N, = 0) on the ensemble. The latter represents a special case (subset) of certain 
symmetry classes in our classification. 



1.1 A few applications of non-hermitean random matrix theory 

As a further motivation to this work, let us give a brief list of physical situations where 
non-hermitean random matrix ensembles play a role. 

An interesting example occurs in condensed matter physics and is due to Hatano and Nelson 
[T5] . A new phenomenon related to a non-hermitean random Hamiltonian is represented 
by the existence of extended eigenstates of the latter - referred to as delocahzation - in one 
and two spatial dimensions. The delocahzation transition can be mapped to the depinning 
of fiux lines from columnar defects in a superconductor at a critical value of the magnetic 
field perpendicular to the defects. 

A problem with similar aspects was considered by Guruswamy, LeClair and Ludwig in [16] , 
where disordered fermions in 2d were studied using complex random potentials. 

In mesoscopic physics, charge transport in quasi-ld quantum wires is described using the 
Landauer formalism. A central object in this approach is the transfer matrix. The trans- 
mission eigenvalues related to the latter determine the main observable, the conductivity 
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of the wire. The distribution of transmission eigenvalues is determined, within the random 
matrix approach, by the solution to a Fokker-Planck type diffusion equation named after 
its authors the DMPK (Dorokhov-Mello-Pereyra-Kumar) equation. In the presence of 
absorption or amplification in the wire, this equation becomes untractable. However, one 
can instead study the reflection properties of the wire in a similar diffusion equation. In 
this case the scattering matrix of the system is not unitary and it is parametrized in terms 
of a non-hermitean matrix [T7]. The same formalism is applicable to multiple scattering 
of classical waves in a random medium, where one is interested in the statistical properties 
of the transmitted and reflected waves, and the phenomenon of localization likewise takes 
place [18]. 

Random matrix ensembles with complex eigenvalues are also relevant in the field of quan- 
tum chaotic scattering, where universal aspects of bound state statistics and resonance 
poles in open systems are described using random matrix techniques. The pioneering work 
by Verbaarschot, Weidenmiiller and Zirnbauer can be found in reference [IS]; see also the 
works by Fyodorov and Sommers, [131 120]. There are a good number of published studies 
in this field, see for example |21j . 

One of the possible applications of random matrix models relevant to QCD involves substi- 
tuting an appropriate random matrix for the Dirac operator (fermion determinant) in the 
euclidean partition function of a gauge theorj0. The integration over random gauge fields 
can then effectively be done by averaging over the matrix ensemble, which can be achieved 
by choosing a simple gaussian distribution of the matrix elements. In the case of QCD, 
the random matrix partition function then describes the static modes of the Goldstone 
degrees of freedom in a finite volume at non-zero quark condensate. This represents a 
scaling region (the epsilon regime or microscopic regime) where the partition function is 
fully determined by global sjTiimetry [23] . 

An important research field in which a better understanding of the tools available to 
calculations in non-hermitean random matrix theory could potentially be of relevance, is 
in the theoretical investigation of finite-density QCD. The notorious sign problem - the fact 
that the euclidean lattice action used in the probability measure for configuration sampling 
in lattice simulations acquires a complex value for finite chemical potential - effectively 
prevents Monte Carlo simulations in this region (see [21] for reviews). In the random matrix 
approach, by including a chemical potential for baryons, the previously anti-hermitean 
Dirac operator acquires a hermitean piece, which renders the Dirac spectrum complex. 
At present, due to the failure of lattice simulations, random matrix theory serves as an 
important tool of investigation of finite-density QCD, and several analytical results have 
been obtained concerning the correlation functions of ensembles with complex eigenvalues 

^The spectrum of this operator is of interest for spontaneous chiral symmetry breaking in QCD because 
of its intimate link to the quark condensate [22] . 
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(see for example [25] and references therein). QCD at finite chemical potential has been 
discussed by several authors [261 HZ] and is a hot topic in random matrix theory as well as 
in lattice simulations. 



1.2 Analytical work on non— hermit ean random matrix theory 

One of the first to do analytical work on non-hermitean random matrix ensembles was 
Ginibre [28], who made a partially successful attempt to determine, in the mid-1960's, 
the eigenvalue distributions and correlation functions for three non-hermitean ensembles 
with real, complex, and real quaternion matrix elements. His work was later completed 
in the orthogonal case by several authors, some very recently, among which Lehmann 
and Sommers, Edelman, Kanzieper and Akemann, Forrester and Nagao, and Borodin and 
Sinclair [29] . 

After Ginibre's work, several new techniques were developed for analytically studying var- 
ious aspects of ensembles with complex eigenvalues. In the paper by Janik, Nowak, Papp 
and Zahed [30] analytical techniques for non-hermitean random matrix ensembles are re- 
viewed. Over the years, several interesting directions of research have been pursued, with 
a peak of activity around 1997. 

Feinberg and Zee developed the method of hermitean reduction (hermitisation) for non- 
hermitean random matrix ensembles ^31j. Green's functions methods were employed by 
Janik, Nowak, Papp and Zahed to obtain the distribution of eigenvalues in the complex 
plane; see for example [32] • The distribution of complex eigenvalues for symplectic en- 
sembles of non-hermitean matrices was also studied by Kolesnikov and Efetov [33], while 
Zabrodin has studied the growth of the domain of support of the complex eigenvalues with 
growing matrix size [34j . 

Fyodorov, Khoruzhenko and Sommers studied the cross-over from hermitean random ma- 
trices of Wigner-Dyson type to non-hermitean random matrices of the type studied by 
Ginibre, using the technique of orthogonal polynomials [35]. Akemann developed the or- 
thogonal polynomial technique for computing correlation functions for strong and weak 
hermiticity in the spirit of Fyodorov, Khoruzhenko and Sommers, discussing also appli- 
cations to the Dirac operator spectrum in QCD with a non-zero chemical potential, see 
for example [26]; similar studies have been performed by Kanzieper [36]; see also the joint 
work on Ginibre's orthogonal ensemble, [37] . New developments in random matrix theory 
(initiated by Kanzieper and further studied by Splittorff, Verbaarschot, and others), also 
relevant for non-hermitean ensembles, include the connection between replica methods and 
Toda lattice hierarchies, see e.g. [55] . 
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Analytical results for non-hermitean random matrix models have been compared to lattice 
data in several publications. We cite a few (fairly randomly selected) in reference [39] . 

Likely, there will be authors that have been overlooked or forgotten in the above summary 
and we extend our apologies to these. 



2 Random matrix theory and symmetric spaces 

It is not in the scope of this article to give an introduction to random matrix theory. There 
are many good introductory texts on the subject (see for example the book by Mehta [lO] . 
the book by Efetov [H] or one of the reviews [9.])- We will be very sketchy here and 
provide a minimum of detail in order to remind the reader of how the symmetric space 
picture arises in the first place. For the interested, the symmetric space connection has 
been discussed in great detail for the traditional hermitean ensembles in the review 
by Caselle and the author. In order not to disrupt the general discussion to review the 
symmetric space approach, we have collected all the relevant material in appendix Rl 

We reproduce here in Table [T]the correspondence between hermitean random matrix ensem- 
bles and symmetric coset spaces of Lie groups that was discussed and thoroughly explained 
in [10\. Let us again clarify the meaning of this table. 

A random matrix ensemble modeling a Hamiltonian Ti. is identified with the coset space 
iP of a semisimple Lie algebra G* expressed as G* = K © iP where K is a symmetric 
subgroup (see appendix |Al). P or iP is the tangent space of a coset space of nonzero 
(positive or negative) curvature. This tangent space is itself a symmetric space of zero 
curvature, and its group of invariance is non-semisimple because P in that case is an 
abelian ideal, [K, P] C P, [P,P] = 0. While all the ensembles are invariant under the 
symmetry transformation h —>■ khk~^ (where k & K, the subgroup spanned by K), the 
invariance of the (gaussian) random matrix partition function in this case is 



The translation invariance is what leads to the non-semisimple structure. In the table we 
have not explicitly written down the coset P, but in the column labeled we have 

listed hermitean random matrix ensembles (Gaussian or Wigner-Dyson G, chiral P 
wave P, and Bogoliubov-DeGennes B) that model Hamiltonian^ with different kinds of 
symmetries. 

^In the chiral ensembles we model instead a Dirac operator - what matters is the larger symmetry 



h khk-^ + h' 
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In the columns labeled "GjK (G)" and "G* jK (G^/G)" we have listed the pairs of posi- 
tively and negatively curved symmetric coset spaces that correspond to P, i.e. e^ and e*^. 
They are random matrix ensembles modeling scattering matrices and transfer matrices, 
respectively, and these ensembles are listed with their traditional names in the columns 
labeled "X"*"" and [C for circular, 5* for S-matrix and T for transfer matrix ensem- 

bles). We can view the scattering matrix of a physical system with Hamiltonian 7i as e*^ 
[2] (at least in a finite range). Therefore the ensemble of the scattering matrix and the 
ensemble of the Hamiltonian of a physical system are listed on the same line in the table. 
Since there is no such simple relationship between the Hamiltonian and the transfer matrix 
of a given physical system (in particular it is not e^) the ensemble of the transfer matrix is 
not on the same line. But it is nevertheless somewhere in the table in the column labeled 

In the columns labeled "mo", "mi" , "m^" we find the multiplicities of the restricted ordi- 
nary, long and short roots of the root lattice pertaining to the symmetric space. This root 
lattice is in many cases different from the one inherited from the complex extension algebra 
G*^. Note also that the total compact Lie algebra G is the sum of the algebra subspace P 
and the symmetric subalgebra K, and the non-compact algebra G* is obtained from it by 
the "Weyl unitary trick": 

G = K©P, G* = K©iP (2.2) 

The algebra K spans the subgroup of invariance of the random matrix ensemble. Inspec- 
tion of the table shows that for the Wigner-Dyson ensembles with /5 = 1, 2, 4 it is SO{N), 
SU{N), and USp{2N) - a well-known result. We note also that the root multiplicities 
rrio, mi, of the ordinary, long and short roots of the restricted root lattice belonging to 
a particular symmetric space, together with the sign of the curvature of the corresponding 
manifold, fully determine the Jacobian of the transformation to eigenvalue space in the 
random matrix ensemble (equation (1A.7P in appendix [X]). The empty spaces in the table 
correspond to random matrix ensembles whose physical applications have not to our knowl- 
edge been discussed in the literature so far, but they are fully defined in the symmetric 
space picture. 

Physical applications of each of the ensembles have been discussed in Part II of |1U] . 
except T4 J 4 that was obtained later, in [12] , as the transfer matrix ensemble of a carbon 
nanotube with symplectic symmetry and an odd number of transmitting channels. In this 
paper, physical observables pertaining to the carbon nanotube - the conductanc^ 5g, the 

group. 

Since the particular carbon nanotubes studied in |12] have a permanently conducting channel due to 
an unpaired transmission eigenvalue, 5g is actually the correction to the conductance due to this channel. 
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localization length ^ describing its exponential decay, and its variance describing universal 
conductance fluctuations - were obtained using the mapping to a symmetric space and 
expressed entirely in terms of group theoretical quantities. Specifically, 



{\n6g) 



2s 



{mi + ms/2) 



7 



mi + ms/2 



Var(ln 5g) 



2 



mi + ms/2 



where s is the dimensionless length of the carbon nanotube expressed in terms of the mean 
free path / of the conduction electrons, and 7 is a constant depending on the total number 
of eigenvalues and on the multiplicities m,,, m^, of the ordinary, long and short roots, 
respectively [T2] . 



3 Symmetries defining random matrix ensembles 
3.1 Summary of results 

In this paper we take a first step towards the identification of random non-hermitean 
matrix ensembles with symmetric spaces. Clearly, the hermitean ensembles exhaust the 
Cartan classification table of symmetric spaces parametrized in terms of real coordinates. 
We will go beyond this classification in listing a number of non-hermitean matrix ensembles 
defined by certain symmetry criteria. The main result of this work is to show that also these 
ensembles define symmetric spaces, but these new spaces will be complex manifolds, in line 
with the fact that the corresponding sets of random matrices have complex eigenvalues. 

We give explicit representations in Tables [2]-[T0] of a large number of random matrix en- 
sembles. The latter are identified with the set P in the tables, whereas the symmetric 
subgroup in each case is spanned by K. The algebra subspace P and the subalgebra K 
satisfy the commutation relations [K, K] C K, [K, P] C P, [P, P] C K and their sum 
spans the total algebra G or G* in the notation of appendix [A] and equation (12. 2p . 
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Table 1: Irreducible symmetric spaces and some of their random matrix theory realiza- 
tions. Some random matrix ensembles with known physical applications are listed in the 
columns labelled X"*", X° and X~ and correspond to symmetric spaces of positive, zero 
and negative curvature, respectively (for references see [lO])- Extending the notation used 
in the applications of chiral random matrices in QCD, where u is the topological winding 
number, we set z/ = p — q. 



Restricted 
root space 


Cartan 
class 


C/K (G) 


C*/K {C^/C) 




mi 








X- 


An-1 
An-1 
An-1 


A 

AI 

All 


SU{N) 

SU{N) 
SO{N) 
SU{2N) 
USpi2N) 


SL(N,C) 
SU{N) 

SL(N,R) 
SO{N) 

SU''{2N) 

USp{2N) 


2 
1 
4 












^2,0,0 
^1,0,0 
^4,0,0 


pO 

^2,0,0 

po 

^1,0,0 
pO 

^4,0,0 


^2,0,0 

Ti,o,o 

^4,0,0 


BCq {p>q) 


AIII 


SU(p+q) 


SU{p,q) 


2 


1 


2u 




X2,l,2u 




Cq ip=q) 


SU{p)xSU{q)xU{l) 


SU{p)xSU{q)xU{l) 


o + 
'-'2,1,0 


^2,1,0 


Bn 


B 


S0{2N + 1) 


SO(2N+l,C) 
SO(2N+l) 


2 





2 




pO 

2,0,2 




Cn 
Cn 


C 
CI 


USp{2N) 

USp{2N) 


Sp(2N,C) 
USp{2N) 
Sp{2N,R) 


2 
1 


2 
1 






B+ 


B° 

-°2,2,0 

B° 

1,1,0 


^2,2,0 
Tl,l,0 


SU(N)xU{l) 


SU{N)xU{l) 


BCq ip>q) 


CII 


USp(2p+2q) 


USp{2p,2q) 


4 


3 






X4,3,4u 




Cq {p=q) 


USp(2p)xUSp(2q) 


USp{2p)xUSp(2q) 




'^4,3,0 


Dn 
Cn 
BCn 


D 

Dlll-e 
DIII-o 


S0{2N) 

SO{4N) 


SO(2N,C} 
SO{2N) 
SO*{4N) 


2 
4 
4 



1 
1 





4 


R + 
-^^2,0,0 

B+ 

4,1,0 


B° 

-^2,0,0 

B° 

4,1,0 

pO 

4,1,4 


^2,0,0 
^4,1,0 
'^4,1,4 


SU{2N)xU{l) 
SO{4N+2) 


SUi2N)xU{l) 
SO*{4N+2) 


SU(2N+l)xU(l) 


SU(2N+l)xU(l) 




Bq (p>q) 


BDI 


SO(p+q) 


SO(p,q) 


1 





V 




Xl,0,iy 




Dq {p=q) 


SO(p)xSO{q) 


SO{p)xSO{q) 




Tl,o,0 
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The ensembles will be defined according to certain physically motivated symmetry criteria 
described in detail below. Among the ensembles listed in Tables [2]-[T0l we find as special 
cases the ten hermitean ensembles of Table [T] corresponding to an even number of eigenval- 
ues (this excludes Cartan classes B and Dlll-odd) and the three Ginibre ensembles. We 
choose for simplicity even-dimensional matrices and we fix the chiral symmetry operator 
to have an equal number of positive and negative eigenvalues. 

Let us now describe the symmetry criteria we have used to obtain the explicit representa- 
tions of the ensembles. 



3.2 The symmetry criteria 

We use as a basis for our classification the symmetry criteria defined in the paper by 
Bernard and LeClair, |14j. The authors of this paper provide a classification of non- 
hermitean random matrices in terms of these symmetries. However, their paper does not 
provide explicit matrix representations for the ensembles of random matrices (they do, how- 
ever, list inequivalent representations of the matrices defining the symmetries). Though 
not very explicit in their papei0, what they classify is cosets of an algebra, but they are 
referred to loosely as "groups" throughout the paper. In this paper we give explicit repre- 
sentations of all these matrix ensembles and show that all the ensembles classified in [H] 
represent algebra subspaces of the type P in equation (]A.9P in the Appendix, reproduced 
below for convenience. The remaining algebra subspace is a symmetric subalgebra K. As 
already mentioned, the two algebra subspaces satisfy the commutation relations 

[K,K]CK, [K,P]CP, [P,P]CK (3.3) 
defining K clS cL symmetric subalgebra. 

We will now, following [1^, define the discrete symmetries we used in defining the cosets 
P. We also keep the notation of this reference and reproduce the definitions used by 
Bernard and LeClair here. It is desirable to limit the number of symmetry classes under 
consideration by some physically motivated conditions. The resulting classification will 
be complete in the sense that all the ensembles defined by the symmetry criteria to be 
defined below will be classified in tables [2] through [101 Following [H] we consider unitary 
transformation matrices p, c, q and k such that a complex matrix h is transformed according 
to 

''On page 3 they do, however, refer to the operator 7i as a "doubled Hamiltonian" , which shows that 
this is what they had in mind. 



9 



h = -php~\ h = ecch^c'\ h = qh^q'^, h = kh*k-^ (3.4) 
pp^ = cc^ = qq'^ = kk'^ = 1 (3.5) 

where ec = ±. Introducing similar signs in the relations for q and k symmetries does not 
lead to new ensembles, because they can be obtained by the substitution h ^ ih [H]. The 
above symmetries will be referred to as P-type, C-type, Q-type and K-type symmetry, 
respectively. 

The physical motivation for these symmetries is as follows. 

The P-symmetry is chiral symmetry {h,p} = with p playing the role of 75. Requiring 
that the p-action be of order two (i.e. that by applying p twice we get back the same 
operator h) implies = 1. Following [H], we will require that all the symmetries are of 
order two. 

The C-type symmetries are the usual anti-unitary symmetries for hermitean operators. 
For ec = + and hermitean h it is equivalent to time reversal symmetry, [h, cK] = where 
c is unitary and K denotes complex conjugation. Since cc* = ±l|f| it follows c^c^^ = ±1. 
For ec = — and hermitean h, the C-type symmetry is equivalent to particle-hole symmetry, 
{h, cK} = (see for example [5| where this symmetry is stated in the form X = iTi = 
— Sa-X^Sa,, = cr^ (8> 1 for the Bogoliubov-De Gennes Hamiltonian Ti). 

The physical requirement that c^c~^ = ±1 can be obtained by requiring that the action of 
c is of order two: 

h = ecc{ecch^ c-^Y c-^ (3.6) 

This implies c~^c^ = ±1, which corresponds to choosing c symmetric or antisymmetric, 
respectively. Note that these two cases correspond in the classical ensembles to presence 
or absence of space-rotation invariance or alternatively, integer or half-odd integer spin. 

'''Time reversal invariance arises from the fact that (in the absence of magnetic fields), if ^(r, i) is a 
solution of the Schrodinger equation, so is 5'*(r, — t). We define the action of the (anti-unitary) time- 
reversal operator T on a state VE* by T'^ — c^E**, where c is a unitary operator. T should reverse the sign 
of spin and total angular momentum. This requirement can be satisfied by the choice c = e*'^'^" for spin 
rotation or c = e*'^'^" for space rotation, with a standard representation of spin or space rotation matrices. 
Now = cc* = e^"'^" = (— )" = ±1. These two cases correspond to integer and half-odd integer spin 
(or presence or absence of space rotation invariance in case oi Jy). 
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The Q-type symmetry reduces to the statement that the ensemble is hermitean if g = 1. 
Therefore we identify all such ensembles with one of the Cartan classes of symmetric spaces. 
Note also that for hermitean ensembles this symmetry becomes trivial, and so it is only 
meaningful for non-hermitean ensembles, where it imposes conditions of (anti-)hermiticity 
on the random matrix or on its sub-blocks. The transformations are unitary, so it follows 
by applying the same transformation twice that q^^q^ = +1 (in this case we get, because 
of unitarity, = 1). 

For hermitean matrices the K-type symmetry is redundant - it is the same as time- 
reversal symmetry - but for non-hermitean ensembles it imposes different kinds of reality 
conditions on the ensembles. From the requirement that the symmetry be of order two we 
obtain kk* = ±1, which implies that k is symmetric or antisymmetric, respectively. 

In [13], if two matrices h were related by a unitary similarity transformation, they were 
considered to be equivalent. We will keep this rule, which is at the heart of the coset space 
structure, but we will report in the comments to the entries in Table [2] through Table [10] 
and in tables [11] and [12] a few equivalent forms of some of the ensembles (all related by 
unitary transformations), for the following reason. As we compared our tables to some 
ensembles used in the literature, we noticed that sometimes a unitarily equivalent form 
of an ensemble was used, which makes its identification in the tables more difficult. We 
therefore included some of the obviously equivalent forms corresponding to symmetry im- 
plementations compatible with the commutativity constraints (see below). Nevertheless, 
the identification of ensembles may not be immediate. For the PQC-symmetric ensembles 
however, the equivalence relations between different ensembles are transferred to an equiv- 
alence relation for the sub-blocks of q and c. Therefore, some of the equivalent ensembles 
can be read off from the corresponding equivalence relations for non-chiral ensembles and 
are not explicitly reported. 

According to the above, we consider h and h' equivalent whenever 

h' = uhv) (3.7) 

with u unitary. This translates into similar equivalence relations for the symmetry imple- 
mentation matrices p, c, q and k defining the ensemble |14j : 

p' = upu\ d = ucu^, q' = uqu\ k' = uku^ (3.8) 

where we have used the fact that uu^ = u^u = 1. Therefore, two ensembles in the ta- 
bles characterized by a subset of {p, c, q, k} and {p', c', q', k'} respectively, are unitarily 
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equivalent if the subsets are related by equation (I3.8p . 



In case of just one symmetry, we consider the unitarity of p, c, g, k together with the 
constraints = 1, cc* = ±1, = 1, kk* = ±1 to obtain all the possible inequivalent 
forms of these matrices. Let us first consider the matrix p. It satisfies pp^ = p^p = p"^ = 1, 
which means that it is diagonalizable with eigenvalues ±1. The solution p = 1 implies 
h = 0. We will, like in [H], for simplicity consider only the case of an equal number of +1 
and —1 eigenvalue^, not least for the purpose of limiting the number of ensembles, and 
we fix p = diag(l, 1,...,— 1). This means that all our chirally invariant ensembles 

will have the usual block-structure, c and k can be symmetric or antisymmetric, while the 
constraint on q excludes an antisymmetric q. The possibilities are listed in Table [2l Note 
that for ensembles 8 and 15 we have chosen a form for k that displays the real quaternion 
structure of the matrix elements (see the specific comment in subsection 14.21 and appendix 

ID. 

If two or more symmetries are present, we will demand, following [Uj, that the symme- 
tries commute. For example, if we have PC symmetry, we demand —p{eccK^ c~^)p~^ = 
^cc{—php~^)^ . For all the pairs of symmetries we then obtain the commutativity con- 
straints of [Hj, we reproduce below: 

c = ±.pcp^ = e^ppcp^ 1 p = ±ikp*k\ 
q = ±pqp^ = epqpqp\ q = ±cq*c^ = ecqCq*c\ (3.9) 
q = ±kq*k\ c = ±.kck^ 

(here ecp, e^, and are signs). The constraints fl3.9p can be used to determine all the non- 
equivalent forms of the symmetry implementation matrices for any subset of symmetries. 
Starting from p, 




(3.10) 



we can then determine the possible forms of c, q and k compatible with eqs. (13. 9p . As an 
example of this procedure, let's consider an ensemble with PC symmetry. For p = p'^ , (13. 9p 
gives either [p, c] = or {p, c} = 0. The former leads to a block-diagonal c and the latter 
to a block-offdiagonal c. Also, c is always either symmetric or antisymmetric. Inspection 
shows that all the antisymmetric forms 

^This choice has no bearing on the symmetric space structure. That is, for an unequal number of 
positive and negative eigenvalues, we obtain similar symmetric spaces. 
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1 
-1 



/ 1 

-1 



1 

-1 oy 



-1 



1 



-1 



1 
-1 



\ 

1 
-10/ 

(3.11) 



are related by unitary transformations. This does not mean that only one of them has to 
be considered in a general ensemble, as the presence of other symmetries may prevent the 
ensemble from being equivalent to another ensemble with the same type of symmetries, 
because also these other symmetries have to transform according to (13.81) with the same 
unitary transformation that relates c to c'. Thus, different equivalent forms for c leave 
open the possibility of the corresponding ensembles being equivalent, but only if the other 
symmetries of the ensembles are transformed into each other by the same transformation. 
There are also three equivalent possible symmetric forms of c 



In the same way we find two equivalence classes for k which are identical to the ones above. 
For q we find the two equivalence classes 



.-1 o. » ■ » n (3.13) 



^ ^ -1 J ' \1 
For PC symmetry, it turns out that the two inequivalent forms 



give rise to the same ensembles, but with opposite sign of in the two cases (table entry 
11). Entries 10 and 11 have unitarily equivalent forms of c, but p does not transform into 
itself under any unitary transformation taking the two different antisymmetric forms of c 
into each other, so they represent two inequivalent ensembles. 

To summarize, for PC symmetry we can have either a block-diagonal or a block-offdiagonal 
c, and in addition = ±c which would give rise to four distinct ensembles, but because 
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the inequivalent forms in the equation above happen to give rise to the same ensembles, for 
PC there are only three different ensembles instead of four, and they are listed in Table [31 

It is also noteworthy that the classified symmetry combinations exhaust all the symmetry 
classes, as having more symmetries of the same type does not lead to anything new. For 
example, two C-type symmetries generate a P-type symmetry if the product of their ec's 
is —1, etc. Also, combinations of the type QK or CK leads to one of the already listed 
classes, since two of the symmetries C, Q, and K combine to give a third type of symmetry. 
We refer the reader to pH] for more discussion on this point. 

Having now explained through an example how ensembles are distinguished, we consider 
this issue clear and we will not go through the details of all the classification. We refer to 
Tables [2] through [10] for a listing of all inequivalent ensembles. 

The integration measures in the random matrix partition functions corresponding to the 
ensembles will not be discussed in detail here - they should be defined by restriction of 
the invariant group measure to the respective cosets. A gaussian invariant measure is 
the simplest choice and is usually sufficient, since choosing a more complicated polynomial 
integration measure often does not change the universal properties of correlation functions. 
We expect an intimate relationship might exist between the non-hermitean Jacobians and 
restricted root lattices, like for hermitean ensembles. We will also not even make any 
attempt here at identifying the new ensembles with coset spaces of groups, so in the last 
column of Tables [2] through [10] we therefore indicate only where to find the ensembles 
known from the classifications of Cartan and Ginibre and leave the other spaces empty, 
indicating that the particular ensemble is of a new, non-hermitean type. 



4 The matrix ensembles 

4.1 General comments on the tables 

In Tables [2] through [10] we list in the first column the types of symmetry defining the 
ensemble of random (non-) hermitean matrices. All the entries marked with an asterisk 
in this column are commented on in subsection 14. 2[ and the lines are numbered for easy 
reference. 

In the second column we list the explicit implementation of the matrices p, c, q and k 
defining the symmetries in accordance with equations (13. 4p and (13.91) . Often two alterna- 
tives are given. Whenever there are two sign alternatives on the same line of this column. 
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the two alternative signs in the following columns correspond to these two alternatives in 
the same order. 



The combinations of symmetries are constrained by the commutativity constraints (13. 9p 
in case of two or more symmetries. In case a C-symmetry is present, the sign of ec defines 
whether it is a time-reversal type symmetry (ec = +) or a particle-hole type sjTumetry 

In some cases there are other implementations of C-symmetry that are possible (and com- 
patible with commutativity constraints) than the ones given in the tables. For example, 
both for ensembles 4 and 10, a possible representation for c (compatible with the commu- 
tativity constraint for p and c for the ensemble with both P and C symmetry) is 



/ 1 
-1 



1 
-1 0/ 



This structure of c (or k) is the one that leads to a representation in terms of real quater- 
nions in the ensembles having this possibility of choice of basis. However, such a c leads 
to a matrix representation in terms of complex quaternions in the ensembles that do not 
have this possibility. Since the latter do not form a field, we prefer the form of c given in 
the table. We will follow this rule throughout the tables whenever there is an ambiguity 
of this kind. 

As we have seen, there is some margin of flexibility and ambiguity in the implementation 
of symmetries due to the possibility of unitarily equivalent representations. We will to a 
large extent follow reference [Hj in the symmetry classification, but in certain cases we 
deviate from it where we find this motivated (in particular, when we want to display the 
real quaternion structure of an ensemble known to possess it, we may choose to display a 
unitarily equivalent ensemble to the one listed in [E]). This does not of course change any 
of its symmetry properties. It should also be mentioned that we have two more ensembles 
with respect to reference [13], ensembles 29 and 30. Probably these were just forgotten 
in [H] - the corresponding combination of signs of ecp, Cpg, Ccg (+ H — ) was included on 
line 2 in the list of equivalence relations at the bottom of p. 5 of the electronic version of 
[H] , but apparently the corresponding symmetry class was not listed. Also, our tables are 
more redundant than the listing in [H| as we include two sign alternatives on almost every 
line, some of which are unitarily equivalent. 
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In the columns labeled "P" and "K" we list the explicit form of random matrix ensem- 
bles with the given symmetries and the symmetric subalgebra defining the subgroup of 
invariance K, 

h khk-^ {k e K) (4.15) 

of the ensembl^. The random matrix ensemble is always identified with the algebra 
subspace P, while the symmetric subalgebra K obtained from eq. (1A.9I) spans the subgroup 
of invariance. Note that we do not distinguish between subspaces P and iP in our tables, 
as they are equivalent, and both of them are reported in the same column. We stress that 
the column "P" contains algebra cosets corresponding to random Hamiltonians H. or iTi 
(or to Dirac operators in case of chiral ensembles) - not scattering or transfer matrices. 
The latter should be obtained by exponentiating P or iP to obtain symmetric coset spaces 
of Lie groups (see section [2] or [10]). 

In case the column corresponding to K is empty, or the corresponding entry is identical 
to the one in the column for P - possibly after using the Weyl unitary trick P ^ iP - it 
means that P is an algebra, i.e. [P, P] G P. This happens whenever the symmetric space 
is a group manifold and not just a coset space, for example for Cartan classes A, C, D (B 
is not included because we limited ourselves to even-dimensional representations). 

We observe that in tables [2] through [10], ensembles corresponding to opposite values of ec 
correspond to the same symmetric subalgebra and are in this sense dual to each other. We 
find that opposite values of ec correspond to random matrix ensembles P with differing 
structure; however such ensembles are considered equivalent even though not related simply 
by the Weyl unitary trick [H], [12] (Ch. 9). 

Just like for the hermitean ensembles, the algebra subspace P is the tangent space of a 
symmetric space. The difference is that these new symmetric spaces are complex mani- 
folds (i.e. when the random matrix ensemble is non-hermitean) , which can be seen as an 
even-dimensional real manifold on the tangent space of which we have defined a complex 
structure. The discussion of these complex symmetric spaces will be deferred to a future 
publication. We know they are symmetric because K and P satisfy the commutation 
relations ( 1A.9I) defining a symmetric space. This has been checked for each ensemble in 
the process of computing the explicit forms of the algebra subspaces from the defining 
symmetry relations. 

If P is a hermitean Cartan class or one of the three Ginibre ensembles, the corresponding 

^As explained before, for symmetric spaces of zero curvature, i.e. algebra subspaces, the subgroup of 
invariance is actually larger and given by (12. ip . 
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class is indicated in the last column labeled "Symmetry class" where "Ginl" means it is 
the Ginibre ensemble with (3 = 1 etc. The Ginibre ensembles have at most one symmetry 
and are all found in the first table. 

In appendix [B] we have collected the relevant properties of quaternions that are needed for 
the quaternion real ensembles. 



Table 2: Explicit representations of random matrix ensembles having no symmetries or only 
one symmetry. See sections [3] and H] for explanations and comments. In the last column of 
tables [2] through [TO] the Cartan classes known from the classification of hermitean random 
matrix ensembles have been indicated, as well as the three Ginibre ensembles. 



Implementation 



K 



Symmetry 
class 



A 



Gin2 



P 



1 
-1 



A 
B 



X 
Y 



c = 1, e. 



± 



A 

A = ±A^ 



X 

X = -x^ 



1 \ _ , 

-1 r^'~ 



A B 
\C ±A^ ) 
B = tB^, C = tC^ 



X Y 

Z -X^ 



Y = Y^ 



q = 1 



A 

A = A^ 



1 
-1 



( -It ^ ) 
A = Al D = 



X 



X Y \ 
Y^ Z ) 

-xl z = -zt 



k = l 



A 

A = A* 



Ginl 



1 
-1 



1 
-1 



A quaternion real 

(At = A) 



Gin4 
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Table 3: Explicit representations of random matrix ensembles with only PC symmetry. 
See sections |3] and H] for explanations and comments. 



Implementation 



K 



P 



1 
-1 



c = 1, ec = ± 



A 
±A^ 



X 



' X 



V 



1 \ 

-1 ; 
/ 1 \ 
-1 

1 
-10/ 

± 



A B\ 
C D 

) 



\ 



( X Y 
Z -X^ 

U V 
W -U'' J 
Y = Y^, Z = Z^, 

v = v^.w = w^ 









1 


;( 







1 







) 





p 



, Cc = ± 



[I oj' 

A = ±A^, B = ±B^ 



X 
-X^ 



or c 



-l 



, Cc = =F 



4.2 Selected comments on specific ensembles 



In this subsection we comment on specific ensembles. The comments are not exhaustive; 
the same type of comments may be valid for other ensembles than the one commented 
upon. Comments are labeled with the number of the row in the tables on which we find 
the ensemble commented upon. We do not give references to literature where ensembles 
are used, because it would be impossible to quote all the relevant papers. 



8 In [13] this ensemble was represented by 



1 
-1 



Here, the chosen basis corresponds to the one in which the matrix elements are real quater- 
nions. 
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Table 4: Explicit representations of random matrix ensembles with only PQ symmetry. 
See sections |3] and H] for explanations and comments. 







Sym. 


Implement at ion 


p 


K 


Symmetry 
class 


12 PQ* 


/I \ 
^=1 -1 ) 

I 1 \ 
^ ~ i ±1 ) 




( i ; ) 

X = -A't, Y = -yt 


AIII 

(for g = 1) 


13 PQ 


I 1 \ 

^={1 I) 


A = A\ B = B^ 


(o -xt) 





11 The same ensembles correspond to two different implementations, 



or equivalently, 

f 1 \ ( l\ 



12, 16 The two ensembles | y ^ ) j ensembles | ^ ^ ) j 

simply related by the Weyl unitary trick. 

14, 16 The two cases on line 14 are unitarily equivalent to each other, and so are the two 
cases on line 16. 



15 Another form for this ensemble is given by the unitarily equivalent form defined by 



P 



1 








-1 



k 



f 1 
-1 

V 



1 
-1 J 



This was the form used in [T3], but it does not make the quaternion structure evident. In 
this case the subspace P is defined by matrices of the form 
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Table 5: Explicit representations of random matrix ensembles with only PK symmetry. 
See sections |3] and H] for explanations and comments. 



Implementation 



K 



Symmetry 
class 



P 
k 



1 








-1 


1 








±1 


1 








-1 



(It) 

A = ±A\ B = ±B* 



( I ' ) 

X = X\ Y = Y* 



P 



( 1 
-1 



1 
-1 



A, B quaternion real 



( i ; ) 

X, Y quaternion real 



P 
k 



1 
-1 
1 
±1 



A 

±A* 



X 
X* 



I A B\ 

-B* A* 

C D 
V -D* C* J 

and the algebra K by 



Y* 

V 



Y 

X* 



\ 

U V 
-V* U* J 



This example shows how different representations can look; however, the above form does 
give a hint about the quaternion real basis when one recalls that a real quaternion has the 
form 
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Table 6: Explicit representations of random matrix ensembles with only QC symmetry 
(continued in Table [7]). See sections [3] and H] for explanations and comments. 



Sym. 



Implementation 



K 



17 QC 



q = l 
c=l, Cc = ± 



A 

A = ±A* = ±A 



T 



X 

X = X* = -X 



T 



18a QC* 



q = l 
/ 1 \ 
-1 



1 
-1 0/ 



A 

A = A- 
A quaternion real 

(At = A) 



X 

x = -x 

X quaternion real 

(Xt = X) 



+ 



18b QC* 



q=l 
' 1 
, ±1 



\ tB* -A* J 
A = A\ B = 



X Y \ 
X = -xt, Y = Y^ 



19 QC* 



c = 1, ec = ± 



A B\ 
[±B^ C ) 
A = ±A* = ±A^, 

B = tB\ 
C = ±C* = ±C^ 



(?t ^) 

X = X* = -x^, 

Y = -Y*, 
z = Z* = -Z^ 




(On the other hand, for example ensemble 28 for ec = + looks like it may have such a 
basis and it doesn't.) 

18a, 18b These two QC-symmetric hermitean ensembles differ by the sign of Cc. However, 
we have chosen different antisymmetric representations for c in the two cases. This is 
possible since g = 1 so the commutativity constraint for QC is trivial and satisfied for both 
choices of c. Thus, the choice 

^=(±1 I) 
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for c is also possible in ensemble 18a; this choice leads to the alternative form for P, 

( f A B ^ 



The ensemble P'_ associated with 



■.a = a\ b = ±b^ 



^=(-11 I) 

is not expressed in terms of real quaternions, but it is related by a unitary transformation 
to the form given in the table for 



/ 1 

-1 



1 
-1 



The ensemble obtained with c 



1 

1 



on line 18b is the one obtained by Altland and 



Zirnbauer for particle-hole symmetric Hamiltonians of Cartan class D in [5j . Inspection of 
its form shows that it indeed has the same form as the matrices of this class on line 17. 



Alternative forms of ensemble 19 are given in Table [TTl 



Lines 20a, 20b describe ensembles P of self-dual and anti-selfdual matrices, respectively. 
These ensembles, that have been listed on two separate lines for clarity, are related by the 
Weyl unitary trick and correspond to the same symmetric subalgebra. 

Unitarily equivalent representations of the ensembles on lines 21a, 21b corresponding to 
the same K, are given in Table [T2j 

On line 23 for ec = ±, P consists of self-dual and anti-selfdual matrices, respectively. 

On line 24 we see that the two Cartan ensembles CI and Dili (i.e., Dlll-even in the 
notation of pilj) are dual to each other and correspond to the same symmetric subalgebrej^ 
The two ensembles are not related by the Weyl trick, and correspond to different symmetry 
conditions for the submatrix A. That these ensembles really correspond to Cartan classes 
CI and Dili can be verified by comparing the form of the representative matrices with 

®This may at first sight appear surprising since they correspond to Dyson index j3 — 1 and /3 = 4, 
respectively. 
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equations (2.17n') and (2.18n') of Chapter 9 in the book by Gilmore [12]. Indeed, the two 
Cartan classes correspond to the cosets 



SO*(2n)/U(n,C) 

and 

Sp(2n,R)/U(n,C), 

respectively, invariant under the same symmetric subgroup. There are other possible im- 
plementations of the matrices q and c. These are: 



1, c= ( A 5 ^c = T 





^ 








V 1 


I) 


, = =F 



/ 1 

^= -1 



^ ~ I -1 J ' ^~ i -1 



Each of these possibilities corresponds to the same two ensembles given in the table, except 
that the last possibility gives rise to ensembles related by the Weyl unitary trick to the 
original ones. 

On line 25 an equivalent implementation of c is 



On line 26 an equivalent implementation of c is 



^01 
-1 

V 





1 



/ 



On lines 29 and 30, the ensembles look very similar. However, they are never the same, 
because the upper and lower signs correspond to different ensembles. Line 29 corresponds 
to a symmetric c and line 30 to an antisymmetric c, even though the commutativity 
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relations between p, q and c are identical for the two entries of the table. They also 
correspond to different symmetric subalgebras with different symmetries of the submatrices 
Y and V. Apparently, these ensembles were forgotten in [14], except for the equivalences 
of their commutativity relations that were mentioned on line 2 in the tabl^ on p. 5. 



5 Conclusions 



Starting from physically motivated symmetries, we have enumerated a number of hermitean 
and non-hermitean random matrix ensembles defining symmetric spaces, giving explicit 
representations of the ensembles and the corresponding symmetry transformations in terms 
of Lie algebras. Thus, we have proved that just like ensembles of hermitean hamiltonian, 
Dirac operator, scattering and transfer matrices of physical systems with given symmetries 
can be classified as symmetric spaces, so also a large number of non-hermitean random 
matrix ensembles define complex symmetric manifolds. We have tried to include in our ta- 
bles all the inequivalent ensembles defined by combinations of certain physically motivated 
symmetries. Our classification is complete in this sense (but we have not proved there are 
no other possible random matrix ensembles defining symmetric manifolds). 

This constitutes only the first step towards physical applications of the theory of symmetric 
spaces to non-hermitean random matrix ensembles, while such applications can already be 
found in the literature for hermitean ensembles, extending all the way to the computation 
of physical measurable quantities expressed group theoretically in terms of the underlying 
symmetric space and its root lattice. The next step should consist in identifying the 
non-hermitean ensembles in Tables [2] - [10] to understand the nature of the corresponding 
geometrical manifolds. We are particularly interested in the operators corresponding to 
free diffusion on such spaces, as these may lead to immediate physical applications. 

In conclusion, the symmetric space picture provides a new approach to non-hermitean ran- 
dom matrix ensembles and puts them into a unified context with the hermitean ensembles. 

^We also believe there is a typing mistake in the last line of the same table which should read 
corresponding to the ensembles on line 27 and line 28. 
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A Symmetric spaces from a random matrix viewpoint 



A random matrix ensemble is defined by a partition function 

dHPp{H) (A.l) 



wliere if is a random N x N matrix, dH is an invariant (Haar) measure on tlie ensemble, 
and PpiH) is a probability measure containing the random matrix potential V{H): 

Pp{H) oc e-^^^*'-^^^) (A.2) 

wliere c is a constant. The partition function depends on a symmetry index /5, the Dyson 
index, that takes one of the values 1, 2 or 4 for real, complex or real quaternion-valued 
matrix elements, respectively. Its value depends on the presence or absence of time reversal 
symmetry and space rotation symmetry (or alternatively, time reversal symmetry and even 
or odd half-integer spin). Pf3{H) is often taken to have a simple gaussian form, because the 
correlation functions are to a large extent independent of the form of the potential V{H) 
in a certain universal scaling limit, the microscopic limit, in which the dependence on the 
average level spacing has been removed. Ensembles are invariant under an automorphism 
of the type 



H UHU-^ (A.3) 

where ?7 is a unitary matrix belonging to a subgroup of invariance. It is evident that this 
invariance property defines a coset structure, which we will show leads to the identification 
of each ensemble with a symmetric coset space. 

For a given U, the random matrix H can be diagonalized: 



H UAU-^ (A.4) 

where A is (block-) diagonal and contains the random matrix eigenvalues {Aj}. Upon per- 
forming the similarity transformation (IA.4p . the random matrix partition function factor- 
izes: 
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Zp^ jdU j \{dKJp{{K})Pp{{\,}) (A.5) 

i 

This means that the part that depends on the eigenvectors in U can be dropped by choosing 
an appropriate normahzation of Zp, and we are left with an integral only over the eigen- 
values. J(3{{\}) is the Jacobian of the transformation ( 1A.4I) . e.g., for a Wigner-Dyson 



ensemble with index (3 it has the form of a product of Vandermonde determinants 



Jp{{>^^})-X{\K-\f (A.6) 

i<j 

It is the Jacobian that gives rise to the interaction between eigenvalues, because it produces 
a logarithmic repulsive pair potential. In general however, the Jacobian does not have this 
simple form, which is valid only for the Wigner-Dyson ensembles of random Hamiltonians. 
It depends on the type of root lattice and on the curvature of the underlying symmetric 
space. The general form of the Jacobian is given by formulas from the theory of symmetric 
spaces |13] 



^°(g) = na6K+(q-«r" 

^~(g)=naeK+sinh""(q-a) (A.7) 

J^{q)=Y{a&R+ sin'"-(q-a) 

These expressions deserve some explanation as they are at the heart of the correspondence 
between symmetric spaces and random matrices. 

In the random matrix partition function, the matrix H belongs to some semisimple Lie 
group G or, in case of a random Hamiltonian or Dirac operator, to the Lie algebra G of 
such a group, while the matrix U in (lA.4p belongs to a subgroup K d G m. both cases. The 
removal of the unitary subgroup K in (lA.Sp corresponds to taking the coset space GjK^ or 
algebra subspace G/K in case H is in the algebra, as the integration manifold. Both G/K 
and G/K are invariant under conjugation by the subgroup K, and the algebra subspace 
actually has a larger, non-semisimple invariance group p kpk^^ +p' which involves also 
translation, i.e. addition of an element from the algebra P = G/K (we assume that the 
subgroup K is connected). 

The integration manifold then is a (symmetric) coset space, whose radial coordinates are 
associated to the random matrix eigenvalues. These radial coordinates are exactly the 



26 



(q-a)'s in (]A.7p . where q are local coordinates identified with the Aj's and a are restricted 
roots. This can be understood by realizing that each point in the symmetric space is 
conjugate to some element from the maximal abelian subalgebra associated to the space. 



A space is said to be symmetric if there is an involutive automorphism a ^ 1 of the 
Lie algebra G onto itself such that = 1 and a splits the algebra G into orthogonal 
eigensubspaces G = K © P corresponding to the eigenvalues ibl: 



a{X) = +X for X e K, a{X) = -X for X e P (A.8) 

such that 

[K,K]CK, [K,P]CP, [P,P]CK (A.9) 

P spans the tangent space of the coset space G/K = e^. Suppose G = K © P is compact; 
then G* = K©iP is non-compact and defines the dual space G*/K = e*^. By enumerating 
all the involutive automorphisms of the compact real form of a complex algebra, defined 
by the following combination of generators Hi in the Cartan subalgebra and raising and 
lowering operators E±a 



(the corresponding real form has a diagonal metric, where as a metric on the algebra we 
use the Killing form), we obtain all the (non-compact) symmetric spaces associated to the 
complex algebra. A symmetric space is characterized by its restricted root system, which 
is obtained with respect to the maximal abelian subalgebra in P. 

As was reviewed in [10], on a Lie algebra G with basis {Xj} a metric tensor gij can be 
defined by the Killing form K{Xi, Xj): 



g,^ = K{X,, X,) = tr(adX,adX,) = Ciq, (A.ll) 

where ad denotes the adjoint representation and C*^ are the structure constants of the 
algebra. A theorem due to Cartan guarantees that the Killing form is non-degenerate 
for a semisimple algebra. We can also define a curvature tensor and sectional curvature 
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(the same as scalar curvature in 2d). It is easy to show that due to the commutation 
relations (]A.9|) . the symmetric manifolds have either positive, negative, or zero curvature. 
On a symmetric space with zero curvature, which corresponds to the algebra subspace P 
itself, the Killing form is degenerate. We can use the normal Euclidean metric in this case, 

gij = ^ij- 



Let us now discuss equation flA.7|) in more detail, repeating some of the above statements 
for clarity. 

[i] The superscript on J. In eq. ( ]A.7p . the superscript on J (0, ±) denotes the curvature 
of the symmetric space. The space whose points are parametrized by the sets of eigenvalues 
{A,,} of random Hamiltonians, having the structure of generators in the coset space G/K of 
a Lie algebra, has zero curvature. These zero curvature spaces actually have a symmetry 
group larger than K, because they are also translationally invariant. Their symmetry 
group is non-semisimple, being the semidirect product of K with the abelian subalgebra of 
translations, corresponding to addition of an element from the algebra subspace P itself. 
A space of random scattering matrices can (at least in a finite range) be thought of as 
the space of matrices defined by ~ e*^ (where Ti. is the Hamiltonian of the physical 
system). This corresponds to a space of positive curvature, while the spaces with negative 
curvature are transfer matrix ensembles. Because there is no simple relationship between 
the Hamiltonian and the transfer matrix of a physical system, the ensemble of the transfer 
matrix is not the negative curvature ensemble corresponding to the symmetric subgroup K 
of the Hamiltonian and the scattering matrix. It corresponds to another symmetric space 
G*'/K'. 

[ii] The product over a G R~^. The product goes over all the positive restricted roots 
of the root lattice belonging to the symmetric space. We remind the reader that the roots 
{a,} of a Lie algebra are defined as the weights of the adjoint representation 



[H„E^] = aiE^ (^ = l,...,r) (A.12) 

where Hi is in the Cartan subalgebra and Ea is a raising or lowering operator. The rank 
of the algebra r is the dimension of the Cartan subalgebra and the roots a = (ai, ...,ar), 
being functionals on the Cartan subalgebra, are r-dimensional vectors in the space dual 
to the same. The restricted roots are similar objects, but to define these one takes Hi in 
the maximal abelian subalgebra (of dimension r' < r) in the algebra subspace P and Ea 
in the algebra G. The restricted root system thus obtained can be of an entirely different 
type from the root lattice inherited from the complex extension algebra. 

[iii] The quantities q ■ a. The quantities q ■ a are exactly the radial coordinates on the 
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symmetric space. In general, if p is a point in the symmetric space P, p' = kpk^^ G P 
(this follows from the commutation relations for K and P). If K is compact, every p & P 
is conjugate with some element h = where H is in the maximal abelian subalgebra of 
P: 

p = khk-^ /i = eS.'?»^^ (A.13) 

where k E K and h is the collective spherical radial coordinate of the point p E P. The 
radial coordinates are obtained as the nonzero diagonal elements of the adjoint represen- 
tation of the diagonal matrix H = Inh in the maximal abelian subalgebra of P, because 
they are the eigenvalues in this representation, and they consist of opposite pairs of the 
form ±q ■ a where a is a restricted root. The local coordinates q = {q^,...,q^ ) on the 
symmetric space are identified with random matrix eigenvalues. The product q ■ a takes 
the form gj ± qj, qi or 2gj for ordinary, short and long roots, respectively. This is exactly 
how the random matrix coordinates appear in the Jacobians of matrix models. 

[iv] The exponents rria. The exponents in (1A.7P are the multiplicities of the restricted 
roots. We remind the reader that in a general root system there are three types of roots: 
long, ordinary and short roots. There is a finite number of possible root systems, because 
the angle between two roots is limited to a multiple of 7r/4 or 7r/6, which follows from the 
general theory of weights and roots. 

The correspondence between the roots and the characteristic quantities of a random matrix 
ensemble also extends to the Dyson and boundary indices f3, a = rrig + mi (and thereby 
also to the parameters defining the orthogonal polynomials, see [lOj): 

mo = /3, mi = P — 1, ms = (A. 14) 

where u is the asymmetry index equal to the number of zero modes of the Dirac operator in 
chiral matrix models. This follows from the identification of the Jacobians of the ensembles 
with flATTl) . 



B Quaternionic ensembles 

Following reference [H] , where a general-purpose introduction to the properties of quater- 
nions and quaternion determinants can be found, we use the representation 
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1, ei = -i(T2, 62 = -i(Ti, 63 = +ias (B.l) 
where cTj are the Pauh matrices, for the basic quaternion units. They Cj's satisfy 

el = el = el = 616263 = —1 (B.2) 
A (complex) quaternion is defined as 



q = qo + e-q (B.3) 

where q= (gi, q2, qs) and go, Q's are complex numbers. A complex quaternion has a dual, 
complex conjugate and hermitean conjugate, defined respectively by 



q = qo-e-q, q* = ql + e-q*, q^ = q*-e-q* (B.4) 



A real quaternion is defined as in ( IB.Sp but with go, qi real numbers. The complex quater- 
nions do not form a field (that is, the inverse of a complex quaternion may not exist), but 
the real quaternions do [H]. Therefore we may consider matrices that can be expressed 
in terms of real quaternion elements. An N x N matrix with real quaternion elements 
can be written also as a 2N x 2N complex matrix, because each quaternion element is a 
2x2 matrix. The real quaternion structure is evident after a certain choice of basis of this 
2N X 2N complex matrix, i.e. after a unitary transformation that transforms it into a real 
quaternion matrix. 

The dual of a quaternion matrix A = [aij], where a is a real quaternion, is defined as 



A = [a,i\ (B.5) 



It follows from definition flB.4|) that a real quaternion matrix satisfies = A. A symplectic 



matrix is defined matrix for which A A = 1. 

We note in Tables [2] through that for the ensembles that can be expressed in terms of 
real quaternions, the antisymmetric form 
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/ 1 
-1 



V 



1 

-1 oy 



-ei 



(B.6) 



of the k- or c-matrix defining the symmetry condition makes this structure evident. This 
can easily be achieved by rewriting the 2x2 complex submatrices of the ensemble as 
complex quaternions and observing that for any quaternion q 



eiq^ ei = q, -eig^*^ei = q* {q^*^ = q^ + e* ■ q*) 



(*) 



(B.7) 



in the C- or K-symmetry relations. 
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Table 7: Explicit representations of random matrix ensembles with only QC symmetry 
(continued from Table E])- See sections [3] and H] for explanations and comments. 



Sym. 



Implementation 



K 



20a QC* 



Q. = 

( ^1 
-1 



1 
-1 



1 

-1 o; 



( A %B\ 
\iB C ) 
A = A, C = C, 
A,B,C quaternion real 



X = -X, z = -z 
X, Y, Z quaternion real 



20b QC* 



q = 

( \ 

-1 



1 
-1 



1 

-1 o; 



( iA_ B \ 
\-B iC j 
A = -A, C = -C, 
A,B,C quaternion real 



X = -X, z = -z 
X, Y, Z quaternion real 



21a QC* 



/ 1 \ 
~ i ±1 j ' ~ + 



A B \ 
\ tB* ±A* J 
A = Al B = ±B^ 



X 



X Y \ 

ty* X* ) 

-X\ Y = tY^ 



21b QC* 



±10'" 



A B \ 
\ ±B* ±A* J 
A = A\ B = tB'^ 



X 



X Y \ 

ty* X* ) 

^_Xt, Y = tY^ 
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Table 8: Explicit representations of random matrix ensembles with PQC symmetry (con- 
tinued in Table [9]). See sections [3] and H] for explanations and comments. 



Sym. 



Implement at ion 



K 



Symmetry 
class 



22 PQC 



P 



1 
-1 



q = l 



1, e, = ± 



A 
±A^ 
A = ±A* 



X = x* = -X^, 
Y = Y* = 



BDI 



23 PQC* 



q=l 
/ 1 \ 
-1 



1 

-1 o; 

ec = ± 



[±A o) 
A quaternion real 



( ■:■ ) 

X = -X, Y = -Y 
X, Y quaternion real 



Cll 

(for ec = +] 



24 PQC* 



P 



/ 1 
" -1 
q = l 

1 \ 

1 j ' ^' 



± 



A' 
±A* 
^A = ±A^' 



X 
\ X* ) 

X = -xt 



CI 
Dill 



25 PQC* 



/ 1 
^=1,0 -1 

/on' 



1 



A 
±A* 
A = A^ 



.Y 

\ X* J 

x = -x^ 



c= 1, ec = ± 
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Table 9: Explicit representations of random matrix ensembles with PQC symmetry (con- 
tinued from Table [H] and continued in Table flOl) . See sections [3] and H] for explanations and 
comments. 



Sym. 



Implementation 



K 



26 PQC* 



P 



1 








-1 





1 


1 





1 
















-1 




± 


1 








-1 





1 


1 






/ 

±c^ tb^ 



A B\ 

±5t C 

J 



( X Y 

Z -X^ 

_yt X* J 

~ 1^ -W )' 
' X Y 

z -x^ 



u 



27 PQC 



P 



(2 J) 

A = A* = ±A^, 
B = B* = ±B^ 



X 
-X^ 

X = x* 



28 PQC 



( l\ 

/ 1\ 

-1 

1 

V-i J 

ec = ± 



( 



A B \ 

TB* ±A* 

C D 
V tD* ±C* 

A = A^, B = tB^, 
C = C\ D = tD^ 



( X Y 

y * X* 

_Xt Y^ 

~ 1^ -zt J ' 
Z={ ^ ^ ^ 

~ \ -r* X* 
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Table 10: Explicit representations of random matrix ensembles with PQC symmetry (con- 
tinued from Table [H]). See sections [3] and H] for explanations and comments. 



Sym. 



Implementation 



K 



29 PQC 



P 



1 



-1 y 
/no \ 

-1 

1 
-1 / 

1 \ 

1 

1 
10/ 



v 



( 

At ±B^ 



A B \ 

TB* ±A* 

J 



I X Y 

U V 
-V* U* J 
X = -Xt, Y = -F^, 
U = -t/t, V = 



ec = ± 



30 PQC 



P 



1 



[0 -1 J 
(1 \ 
-1 

1 

V 0-1/ 

^01 \ 
1 

1 

V -10/ 



( 

V-5t ±A^ 



A B \ 

±B* ±A* 

J 



( X Y 

U V 
V* U* J 
X = -xt, Y = Y^ 
U = -f/t, v = v^ 



ec = ± 
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Table 11: Unitarily equivalent representations to the ensemble on line 19. The signs refer 
to the sign alternatives on line 19 and a prime denotes yet an alternative form (e.g. line 
19' is equivalent to the two ensembles on line 19; hne 19+ and 19'+ are equivalent to the 
upper sign alternative on line 19; etc.). 



Sym. 



Implementation 



K 



19' QC 






19+ QC 






19- QC 




c = 1, ec 





19'+ QC 



1 

1 



1 

1 



, ec = + 





19'- QC 




1 

1 
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Table 12: Unitarily equivalent representations to the ensembles on lines 21a and 21b. 
Ensemble 21a+ is equivalent to the upper sign alternative in ensemble 21a etc. 



Sym. 



Implement at ion 



K 



21a+ QC 



1 

1 



c — 



1 
-1 



, Cc = + 



A B\ 

\ -B* A* J 



X 



X Y 



-rt 



21b+ QC 



1 
\ 



1 
-1 



(A B \ 

\ B* -A* J 
A = -A'^, B^B^ 



X Y\ 
\ -Y* X* J 

X = -x'^, y = yt 



21a- QC 



_ / 1 
" 1^ 1 
1 ^ 



c = 



-1 



, Cc = + 



A B \ 
[c J 

A = A* 
B = -B* = -B^ 
C = -C* = -C^ 



x = x* 

Y = Y* = Y^ 

Z = -Z* = Z'^ 



21b- QC 



_ / 1 
~ 1 
\\ 



-1 



, ec = - 



I A B 
\C -A^ ) 
A^-A* 

B = B* = B^ 
C = C* = C^ 



[z 

X^X* 

Y = Y* = y^ 

Z = -Z* = z^ 
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